Using CMBR analysis tools for flow anisotropics in relativistic heavy-ion collisions 
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Recently we have shown that there are crucial similarities in the physics of cosmic microwave 
background radiation (CMBR) anisotropies and the flow anisotropics in relativistic heavy-ion col- 
lision experiments (RHICE). We also argued that, following CMBR anisotropy analysis, a plot of 
root-mean square values of the flow coefficients, calculated in a lab fixed frame for RHICE, can yield 
important information about the nature of initial state anisotropies and their evolution. Here we 
demonstrate the strength of this technique by showing that elliptic flow for non-central collisions 
can be directly determined from such a plot without any need for the determination of event-plane. 



I. INTRODUCTION 



Relativistic heavy-ion collision experiments (RHICE) are often termed as Little Bangs in analogy of the Big Bang 
representing the initial stage of the Universe. Indeed, there are tempting similarities between the early universe and 
these experiments. For example, it is often mentioned that the surface of last scattering for the cosmic microwave 
background radiation (CMBR) is similar to the freezeout surface in RHICE in the sense that one has to learn about 
I ■ the early stages of the system from hadrons coming out from the freezeout surface, just as CMBR encodes information 
about the early universe. However, these correspondences and analogies have been primarily made from a motivational 
view point. It has recently been shown by us [1] that the physics of RHICE may have deeper connections, not just 
with the above mentioned aspects of CMBR, but even to the most celebrated aspect of inflationary physics of the 
universe namely the presence of superhorizon fluctuations. 

Following the techniques used for CMBR anisotropy analysis, it was proposed in ref. iJj] that instead of focusing on 
. the average values of the flow coefficients t;„, one should calculate root-mean square values of the flow coefficients 
Q.^ ' v™^. Further, these calculations should be performed in a lab fixed frame, which eliminates the difficulties associated 
^vq with determination of event plane for conventional elliptic flow analysis for non-central collisions. To distinguish from 
the conventional flow coefficients u„ which are defined with respect to the event plane, we will denote, in this paper, 
' , the flow coefficients defined with respect to the lab fixed frame as Vn- The root mean square value is then denoted as 
■ v™^. Plots of (denoted as v™^ in [l|) using model estimates based on initial energy density anisotropies from 
QQ ' HIJING [2] were given there for large range of values of n (with n ranging from 1 to 30) for central collisions and it 
r. was argued that such plots can yield important information about the nature of initial state anisotropies and their 
^ ' evolution for RHICE. 

• ^ , In this paper we demonstrate the strength of this analysis technique by showing that a plot of values of 5™*' vs. 
' n can be used for directly probing various flow coefficients, in particular, the elliptic flow for non-central collisions 
H I 3, without any need for the determination of event-plane. For non-central collisions, a plot of u""" (using model 

- - < estimates from the spatial distribution of initial parton energy density), shows a prominent peak at n = 2, as we will 
see below. Further, the height of the peak increases with increasing value of the impact parameter b, thus allowing 
direct determination of elliptic flow (in terms of w™*). This can be very useful for analyzing huge wealth of data as 
the determination of event plane for each event becomes redundant. 

The correspondence between the CMBR physics and RHICE explored at a deeper level as discussed in [1] may lead 
to new set of interesting phenomena and new techniques for RHICE. We mention that such a connection between 
physics of RHICE and that of inflationary universe was never anticipated earlier, and indeed, at first sight, it looks 
surprising that a concept like superhorizon fluctuation which arises from highly non-trivial, superluminal expansion 
phase of the very early universe could have any relevance for relativistic heavy-ion collision experiments in laboratory. 
Such superhorizon fluctuations in RHICE originate from the fact that in the center of mass frame the thermalization 
happens rather quickly, within about 1 fm. This time is too short for any fluctuations (inhomogeneities) in the 
transverse direction to disappear, by processes of homogenization, which have wavelengths larger than 1 fm. Initial 
parton energy density distribution from HIJING show that transverse fluctuations with wavelengths significantly 
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larger than 1 fm are necessarily present at the time 1 fm even in central collisions. These arise from localization of 
partons inside initial nucleons, as well as from the fluctuations in nucleon coordinates. Analyzing the development of 
flow from such initial large wavelength fluctuations, we had argued that the anisotropics in the final particle momenta 
may show characteristic features of acoustic oscillations of subhorizon modes as well as suppression of modes which 
remain superhorizon at the freezeout stage. We will show below that such features do not mask the peak at n = 2 in 
non-central collisions. 

The paper is organized in the following manner. In section II, we recall the adoption of CMBR anisotropy analysis 
method for the case of RHICE. Here we also show that calculations of ■DJ^™* in the lab fixed frame can directly yield 
information about various fiow coefficients defined with respect to the event plane. In particular, the conventional 
elliptic fiow can be determined in this manner. Section III presents the results of the application of these analysis 
methods for the determination of z;™* for the case of non-central collisions. In section IV we discuss the px depen- 
dence, and in section V we discuss the issue of anisotropic detector acceptance. Section VI presents conclusions and 
discussions. 



II. USING CMBR ANALYSIS TOOLS FOR RHICE 



Traditional analysis of elliptic flow in RHICE [s*] is based on measuring the average value of the 2nd Fourier 
coefficient of particle transverse momentum anisotropy. The average value is determined by carrying out the Fourier 
expansion in the event-plane reference frame. Even for the so-called 'central collisions' one determines event plane for 
very small impact parameter collisions and then carries out the traditional analysis of determination of average fiow 
Q. However, determination of event plane is highly non-trivial making determination of higher flow coefficients very 
difficult. Indeed, flow coefficients have been determined only upto ve- 

It is here that the use of techniques of CMBR analysis can make a significant difference. For CMBR, the temperature 
anisotropics are analyzed using spherical harmonics, as appropriate for the surface of 2-sphere (the CMBR sky) 0]. 
The coefficients of the expansion (denoted as aim corresponding to the spherical harmonic Y/m) are degenerate in the 
argument m and when averaged over different values of m yield zero averages due to isotropy of the universe (when 
suitably corrected for local velocities etc.). What one plots is the variance of a/,„ denoted by C/ and this leads to the 
celebrated power spectrum of CMBR anisotropics ^ . 

We propose using the same technique jlj for analyzing particle momentum anisotropics, using lab fixed frame, in 
RHICE to probe the generation and evolution of flow. For RHICE, focusing on central rapidity region, one will 
be analyzing momentum anisotropics on a circle, requiring use of the Fourier coefficients, which we denote as Vn to 
distinguish from the conventional fiow coefficients w„ which are defined with respect to the event plane. With a fixed 
lab frame, the event average values of these VnS will all be zero due to rotational symmetry. We then propose to use 
the variance of Vn, i.e. v™^ in analogy with C; for CMBR. 

We first establish the relation between calculated in the lab fixed frame and the conventional flow coefficients 
Vn which are defined with respect to the event plane as follows 



E 



(1) 



Where the angle (f>' is measured with respect to the event plane. Assuming that the event plane is oriented by angle 
tp with respect to the lab fixed frame, we can write the above equation in the lab fixed frame as. 



p^4^)= J2 (2) 

71— — oo 

where now the angle (f) is measured with respect to the lab fixed frame. The subscript for p^{4>) remind us that the 
above Fourier series expansion uses the flow coefficients defined with respect to the event plane which is oriented at 
an angle 4' in the lab frame. Without the knowledge of ip, we will calculate the lab fixed flow coefficients Vn{ip) as 



VnW = ^ r p^{<f)e-"''^d<l> (3) 

Here the argument ■0 in Vn{ip) is used to remind that these quantities are calculated in the lab fixed frame, and are 
different from the values of conventional flow coefficients u„ defined in Eq.(l). The event average value of these 
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flow coefflcients Vn, is obtained by averaging {;„ over a large number of events (of similar type, e.g. centrality etc.). 
Vn is determined by the following expression 

^^^f_ P^(<^)e-™*d0)dV (4) 

Integration over reflects the random variation of the orientation of the event plane in the lab fixed frame from 
event to event. By changing the order of integration for d(j) and d^, and using Eq.(2), we immediately see that the 
Vn — Q for all n > as mentioned above. 

We next calculate the event averaged values of as 



27r 

Straightforward calculation using Eqs.(2)-(5) gives 



{vT'f^vl^^l d#„(VK(V) (5) 



vT' = \yn\ (6) 

This is an important result. It shows that various conventional flow coefficients u„, defined with respect to the event 
plane, can be directly calculated in the lab fixed frame in terms of "D™*. However, one has to be careful in interpreting 
the above equation. The above equation has been derived for the event averaged values. Thus the right hand side 
\vn\'^ also represents the event averaged value of the conventional flow coefficients. Even if we assume that each event 
is almost identical in terms of centrality selection etc. (which it is not), still there will be random fluctuations in the 
fluid especially due to initial state fluctuations If the contributions of such fluctuations becomes dominant then 
the right hand side of the above equation will not have a neat interpretation of relating to defined with respect 
to the event plane. 

Fortunately this problem does not exist for the elliptic flow V2 ■ Below we will see, with model calculations of v™'^ 
for non-central events from HIJING, that the contribution arising from the elliptic shape of the produced partons 
in the overlap region is very large compared to the contribution from the random, initial state fluctuations. This 
will be seen as a prominent peak in the plot of i;™* vs. n at n = 2. Thus, as far as determination of the elliptic 
flow is concerned, Eq.(6) shows that calculation of {12""* in the lab fixed frame can be used to determine the value of 
U2I. Neglecting the contributions of random fluctuations, symmetry of the elliptical shape of the overlap region for 
non-central collisions will then imply that \v2\ is the same as V2 conventionally defined as the coefficient of the cos2(j) 
term. 

For other values of n the Eq.(6) has to be suitably interpreted with proper account of event by event random 
fluctuations. As we have emphasized, values of -D™'* determined by averaging over events should be directly used to 
probe the statistical properties of fluctuations and anisotropics of the initial plasma region (as is done for CMBR). 
Hydrodynamical simulations, with proper incorporation of such fluctuations can then be used to predict the values of 
-rms •^yjjose comparison with data can then be used to constrain/determine various physical inputs such as equation 
of state, viscosity etc. 

We now describe our method for calculating w™* using HIJING [6|. We first recall from ref. ^ the estimates of 
initial spatial anisotropics. We start with the initial transverse energy density distribution for Au-Au collision at 200 
GeV/A center of mass energy from HIJING. For parton positions we use random locations inside the parent nucleon, 
(similarly, for partons produced from the string systems, random position are used along the line joining the two 
corresponding nucleons). The transverse energy density at a given transverse position x, at proper time r = Teq, is 
taken as 0, 



i 

where denotes the initial transverse coordinates of the ith parton (determined using the coordinates of the 
parent nucleon in HIJING as discussed above), El^ is its transverse energy, ptr the transverse momentum, and i? 
is its transverse velocity. For the rapidity window we take A(j/') — 1 centered at ?/ = 0, 0. The sum over i 
includes all partons in a small transverse area element Ayl(~ 0.5 fm^) at position x. We have included a factor 
F{Teq,ptr) = 1/(1 + ^ / {PtrTeq)'^) to account for the probability of formation of partons with zero rapidity 
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We assume that the hydrodynamic description becomes apphcable by r = Teq, which we take to be 1 fm and 
calculate the anisotropics in the fluctuations in the spatial extent i?(0) at this stage, where RUj)) represents etr 
weighted average of the transverse radial coordinate in the angular bin at azimuthal coordinate (/) [i| . As emphasized 
above, angle (f) is taken in a lab fixed coordinate frame. We divide the region in 50 - 100 bins of azimuthal angle (/>, 
and calculate the Fourier coefficients of the anisotropics in SR/R = {R{(t>) — R)/R where R is the angular average of 
R{(j)). Note that in this way we are representing all fluctuations essentially in terms of fluctuations in the boundary 
of the initial region. We use Fn to denote Fourier coefficients for these spatial anisotropics, and use Vn to denote nth 
Fourier coefhcient of expected momentum anisotropy in 5p/p defined in the lab frame. Here 5p represents fluctuation 
in the momentum p of the final particles from the average momentum, in a given azimuthal angle bin. 

We should clarify that in the conventional analysis, V2 is directly related to the elliptical shape of the fireball, and 
normally one does not call it a fluctuation. However, in our language, the elliptic shape of the fireball itself is taken to 
represent a fluctuation (of large wavelength) from isotropic case. As our analysis is carried out in the Lab flxed frame, 
all shape fluctuations (including the elliptical ones) are treated uniformly by calculating different Fourier coefficients 

Fn. 

Again we emphasize the most important difference between the conventional discussions of the elliptic flow and our 
analysis; here one does not try to determine any special reaction plane on event-by-event basis. A flxed coordinate 
system is used for calculating azimuthal anisotropics. This is why, as discussed in Sect.H and as is seen in the plots in 
ref.Q, averages of i^„s (and hence of Vns) vanish when large number of events are included in the analysis. However, 
the root mean square values of -F^s, and hence of v^s, are non-zero in general and contain non-trivial information. In 
fact, it is the same as the standard deviation for the distribution of F„s since the average value of F„s is zero. This 
is what is exactly done for the CMBR case also 

One important difference from CMBR, which may be of crucial importance in flnding any non-trivial features in 
the plot of ■OJ^™* for RHICE, is the fact that for the CMBR case, for each I mode of the spherical harmonic, there 
are only 2^-1-1 independent measurements available, as there is only one CMBR sky to observe. For small I values 
this becomes dominant source of uncertainty, leading to the accuracy limited by the so called cosmic variance 
In contrast, for RHICE, each nucleus- nucleus collision (with same parameters like collision energy, centrality etc.) 
provides a new sample event. Also, in the absence of any special reflection symmetry here (which was present in the 
traditional elliptic flow analysis) all flow coefficients give non-zero contributions to w™", and the "sin" terms give 
same values as the "cos" terms. In the plots of ■OJ^™* below, we show the sum of these two contributions, i.e. square 
root of the sum of the squares of the "sm" term and the "cos" term 

III. RESULTS FOR NON-CENTRAL COLLISIONS 

We have generated events using HIJING and we present sample results for Au-Au collision at 200 GeV/A center 
of mass energy. In all the plots, the averages are taken over 10000 events, and the root mean square values z;™* of 
the flow Fourier coefficients are obtained from spatial -F^s simply by using proportionality factor of 0.2. The relation 
between the Fourier coefficients of the spatial anisotropy and resulting momentum anisotropy in our model can only be 
obtained using a full hydrodynamical simulation, with proper accounts of any surface tension, as well as factors such 
as horizon crossing etc. to properly account for the physics discussed here. In the absence of such a simulation, we 
make a strong assumption here that all Fourier coefficients for momentum anisotropy are related to the corresponding 
coefficients for spatial anisotropy by roughly the same proportionality factor, which we take to be 0.2 for definiteness. 
(As mentioned in ref.[l[, this choice also gave reasonably good agreement with the results for -D™* in the literature 
for (almost) central events Q.) Although overall shape of the plot of will crucially depend on the values of 
these proportionality constants, and hence may change completely from what is shown here, the qualitative features 
such as the presence of peaks may remain unaffected. (Though, the peak positions may shift depending on these 
proportionality constants.) F„s are calculated directly from the parton energy distribution at i = 1 fm for each event, 
within unit central rapidity window. 

Fig.l presents plots of v™'' for different values of impact parameters. Solid, dashed, and dotted plots correspond 
to impact parameter 6 = 0, 5 fm, and 8 fm, respectively with the spread in b taken to be 1 fm for non-zero b cases. 
The solid plot for central collisions shows monotonically decreasing values of v™^ . We emphasize that even here the 
plot has non-trivial structure compared to the case when partons are randomly distributed inside the nuclear volume 
with uniform probability, which yields a flat plot for corresponding {j™*. (Though, the rise oib — Q plot for small n 
does not appear prominent here due to larger vertical scale of the plot.) Errors in these plots are very small. 

Important thing to note is the prominent peak in non-central collisions for 6 = 5, and 8 fm at n = 2. The peak 
height is signiflcantly larger for 6 = 8 fm. Important thing to realize is that this feature (which within our model 
simply represents elliptical distribution of partons in position space at the initial time) is present even when the 
v™^ have been calculated in a lab fixed reference frame. Average values of Vn for all values of impact parameters. 
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FIG. 1: These plots show smooth joining of the values of C^'"" calculated from Fourier coefficients F-nS of the initial spatial 
anisotropy as discussed in text. Plots are obtained using 10000 events from IflJING. Solid, dashed, and dotted plots correspond 
to impact parameter & = 0, 5 fm, and 8 fm, respectively. Note the prominent peak at n = 2 for non-zero b cases. 



calculated for large number of events, continue to be zero as shown in ref. [l| for the central collision case. This is 
expected from rotational symmetry in lab fixed frame. Larger overall values of {i™*' for larger values of b are simply 
due to smaller number of initial partons leading to larger fluctuations. We mention that the solid plot in Fig.l is 
given exactly for 6 = to show the difference between zero and non-zero b cases. However, with experimental data, 6 
= will inevitably include events with small but non-zero values of b. Thus even for b = case one should expect to 
observe a (small) peak at n ~ 2. 

We have also calculated the root mean square values of the Fourier coefficients of the momentum anisotropy directly 
using the momenta of final particles from HIJING. These plots are all similar to each other for all values of the impact 
parameter b and hence similar to the central case 6 = (apart from larger statistical fluctuations for larger values 
of b due to smaller number of partons as discussed above for Fig.l.) This is expected because, despite initial spatial 
anisotropy of parton production for non-central collisions, momentum distribution remains isotropic in the absence 
of any hydrodynamic evolution. We do not show these momentum plots here because such plots have been given in 
ref.[l| for the central collision case. 

Fig.l shows the strength of this technique of plotting root-mean-square values of flow coefficients in lab fixed 
coordinate frame as compared to the traditional analysis of elliptic flow. This technique is simple to implement (apart 
from the issue of non-flow correlations 0, Q , which one has to resolve in any case) , allowing larger statistics to be 
generated. What one needs is predictions from hydrodynamics about u""*^ calculated in lab fixed frame and then 
compare with data analyzed similarly. As the plots in Fig.l show, the values of v^"^^ even for very large values of n 
continue to be significantly non-zero and with the possibility of sufficient statistics, one may be able to determine any 
important features in these plots. 

In [l*] we had argued that certain important aspects of inflationary density fluctuations, such as suppression of 
superhorizon modes and acoustic oscillations may also be present in these plots of -D™* in RHICE. Basic ideas can 
be stated as follows. Acoustic peaks in CMBR primarily result from the coherence and acoustic oscillations of the 
inflationary density fluctuations with coherence originating from the superhorizon nature of fluctuations, essentially 
freezing the fluctuations on superhorizon scales. This should be reasonably true for RHICE as the transverse velocity 
to begin with is expected to be zero and becomes non-zero only due to pressure gradients (which docs not become fully 
effective on superhorizon scales). The oscillatory behavior for the fluctuations in the universe results from attractive 
forces of gravity and counter balancing forces from radiation pressure (with the coupling of baryons to the radiation). 
The oscillatory behavior of the fluctuations in RHICE is expected due to flow development with non-zero pressure 
gradients, at least for sufRciently small wavelength modes. These arguments suggest that some features like acoustic 
peaks may be present in the plots of -D""" for RHICE also. Further, modes with wavelengths larger than the acoustic 
horizon size H^^ at freezeout should be naturally suppressed as pressure gradients do not become fully effective on 
those scales to generate full flow anisotropy before freezeout occurs. This leads to the following suppression factor 



{"^n) observed — ^ ('^n)maa; (S) 

where A ^ 2TrR-f^ /n, (n > 1), is the measure of wavelength of the anisotropy corresponding to the nth Fourier 
coefficient. Here R^"^ represents the transverse radius at the freezeout time Tf^- Using the rough estimate of the 
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FIG. 2: These plots are obtained from the values of v™" of the plots in Fig.l by including (a) the suppression factor given in 
Eq.(8) for superhorizon modes, and (b) including acoustic oscillations as well. The plots here only show the modeling of the 
type of suppression factor and oscillations discussed in the text, which are superimposed on the plots in Fig.l. 



rate of change of the transverse velocity to be about 0.1 fni^^ at the early stages at these energies Q, we can 
estimate R^'' ~ ^ + 0.05(r/r — Tg,)^ = ^(1 + Here R = R{Teq) = Cs(r/r — Tgg). The largest wavelength 

Amaa; of Spatial anisotropy which will have chance to develop to its maximum hydrodynamic value is, therefore, 
\max — 211^^ = 2cs{Tfr — Teg) = 2R{Tf.q) . Thls glvcs US the Corresponding minimum value Umin of n below which 
flow coeflicients should show suppression due to being superhorizon. 



The possibilities of the superhorizon suppression and acoustic oscillations lead to new peak structures in the plot 
of v™^ as was shown in [l[ for central collisions. It is therefore important to know whether the peak in ■D™'* due 
to elliptic flow for non-central collisions (as in Fig.l) survives if these new features are also present, especially with 
the suppression of superhorizon modes at smaller n. For this purpose, we show in Fig. 2, plots (for different values of 
impact parameter b as in Fig.l), when superhorizon suppression is included and when acoustic oscillations are also 
included in modeling the values of -0™*. It is clearly seen that the peak dX n — 2 corresponding to the elliptic flow is 
present in all the plots and remains most prominent for 6 = 8 fm. 

Here we mention that recently Sorensen has proposed a novel explanation of azimuthal correlations observed 
at the Relativistic Heavy-ion Collider (RHIC) at BNL in terms of the suppression of the superhorizon fluctuations 
discussed in [l[. One of the plots (Fig. 4) in [Tol| (where suitable subtraction of V2 has been made) shows possible 
suppression of values of the variance for small n, as predicted in [I,]- 

We mention that we calculate v™^ by direct calculation of variances of the distributions of (actually F„s) in 
the laboratory fixed frame. Equivalently one can calculate 'D™'* by calculating the two-particle azimuthal correlation 
functions, as is done for the elliptic flow calculations. To extract elliptic flow coefficient from two-particle azimuthal 
correlations one has to separate non-fiow contributions 0,3 • Two particle azimuthal correlations, which are exper- 
imentally measured, contain contributions from non-flow effects such as jets, resonance decays, HBT correlations, 
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final state interactions etc. Various methods have been discussed to separate out the non-flow contributions to the 
azimuthal correlations 0, Q • Our estimates of v™^ will also contain such non-flow contributions which have to be 
properly accounted for. 



IV. pT DEPENDENCE 

In the conventional analysis of elliptic flow, V2 has a non-trivial dependence on pr of hadrons V2 increases with 
Pt, saturates, and then decreases. Initial increase of V2 with pT is easily determined [l^ [isl. 1 14| . Large pt behavior 
of flow is governed by viscous effects, parton energy loss, as well as by parton coalescence The same physics will 
be applicable in our model as well. For simplicity, we will follow the approach in ref. [12i | where the linear increase of 
the elliptic flow V2 with pT is determined to be 



V2 = 7^{PT - vmr) (10) 



where a characterizes the magnitude of elliptic flow, v is the flow velocity (see below) and niT — \Jp\ -V rri^ is the 
transverse mass. We mention that this linear dependence of V2 on px is expected for the intermediate range of pt- 
The behavior for small pr may be more non-trivial [3]. However, recall from Eq.(6) (and discussion in Sect. II) that 
in the context of our model, -y™* is simply given by the conventional V2-, hence whatever be the dependence of V2 on 
Pt, should show the same dependence. Still, it is useful to directly check the expected pt dependence of w™^ 
in our model by adopting the analysis of V2 in ref. [l3| for our case. This is for two reasons. First, the relevance of 
random fluctuations needs to be assessed at different stages of the calculations to see the effects of event averaging 
and to see whether similar statements can be made about other flow coefficients. (For example, see the discussion 
following Eq.(12) below about the form of the fluid 4-velocity.) Secondly, the analysis for V2 in ref. [l3| is carried out 
to linear order in a. Since we calculate the root mean square values, we need to keep terms of order c? . Though, as 
we will see below, the term will turn out to be only relevant for W4 and will not affect the analysis of w™'*. 

We start by noting that the expression for V2 in Eq.(lO) corresponds to the conventional analysis with the determi- 
nation of event plane. In our approach, the analysis is carried out in the lab fixed frame. For each non-central event 
the magnitude of the flow anisotropy will still be characterized by (say) the above equation. However, the maximum 
flow direction will now vary randomly from one event to another because the event plane is left undetermined. Thus, 
it looks reasonable to expect that {i™* should show the same dependence on pT as V2 in Eq.(lO), as is indeed shown 
by Eq.(6). We will check this directly by adopting the analysis of V2 in ref. [l^] for our case for the estimation of 
~rms ^j^g fixed frame (and keeping terms of order a^). Let us start with the transverse momentum distribution 
of emitted particles [l3] 



PW K exp{ ) (11) 



PTdpTdpzd<j) T 

Here, u{(f)) is the space component of the fiuid fiow 4-velocity in (f) direction, and vP^cj)) is the time component. 
For non-central collisions, the anisotropic fiow can be parameterized as 



u{(j)) = u + 2acos2{<j) - ip) (12) 

Here u is the angular average of maximum fiuid 4-velocity, and a > characterizes the magnitude of elliptic flow. 
Note, that we have introduced here the angle ip characterizing the orientation of the event plane in the lab frame so 
that the angle (p is measured in the lab frame. Also, it is important to realize that here we are only characterizing 
elliptic flow anisotropy of the fluid, without worrying about the presence of general fluctuations. The justification for 
this is that, as can be seen from Fig.l, for non-central collisions in our method it is only ■D™'* which has a significantly 
larger value than the values of ■u™'' for other n. Thus, the presence of general fluctuations (as accounted for in Fig.l) 
will play an important role in determination of these other values of ■OJ^™^ even for non-central collisions. In contrast, 
for -D""" the elliptic anisotropy resulting from non-central collision seems to play dominant role over the random 
fluctuations. Hence the parametrization of Eq.(12) may capture important aspects of the dependence of w™'' on pt 
etc. 

We are interested in calculating the root mean square value (instead of V2 as in ref. (l^l. Thus, we expand 
u^{(j)) = y^u{(f>)^ + 1 to second order in a, and get (with v = u/v?), 
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= M° + 2avcos2{4> + -^-pT3Cos22((/) - V-) (13) 

(u ) 

Using this equation, and Eq.(ll), we get particle transverse momentum distribution as 

^ All + 2«ca.2(0 - ^^(PZ^pil + + ,,,4(0 - ^))(^^^1^S^ - ^)] (14) 

where A is independent of ?/'• We see that this equation is in the same form as Eq.(2). It is important to note 
here that the term here only comes for V4 (requiring higher powers of a for the calculation of w™'*) and V2 remains 
linear in a. 

Following the discussion in Sect. II following Eq.(2), and using Eq.(6), we immediately see that 

^^iPT-vmr) (15) 

Which is the same dependence on px as in Eq.(lO) for the conventional analysis of V2 using event plane determina- 
tion. 

We again emphasize that in our picture, the randomness in the orientation of the event plane w.r.t. the lab fixed 
frame is easy to model by introducing angle Tp. The other intrinsic sources of fluctuations will also contribute to -5™^ 
in our method, but as shown in Fig.l, their contributions appear to remain small (as seen by the plot for 6 = 0). 
(Note again, overall larger scale of plots of -D™* for larger values of b is primarily due to smaller number of particles 
leading to larger statistical fluctuations.) As we had discussed earlier, v™'' for n ^ 2 may be dominated by these 
intrinsic sources of fluctuations. If the contribution of such fluctuations remains sub dominant, then our results show 
that direct determination of V2, in particular its pr dependence ([T^) is possible by measuring -O™* in lab fixed frame 
completely avoiding the determination of event plane. 



V. DETECTOR ACCEPTANCE 



One of the main advantages of using the techniques proposed here is that it makes it much easier to extract the 
information about flow without any need for event plane determination. Thus large amount of data can be analyzed 
improving accuracy. It then becomes important to address the issue of being able to incorporate data with varying 
azimuthal acceptance, as has been emphasized for the conventional flow analysis. Several methods have been proposed 
for accounting for azimuthal dependence of detector acceptance in the conventional flow analysis [1, [lH| . These can be 
suitably adopted to the method proposed here. As mentioned above, due to fixed lab frame being used for analysis, 
the average values of all flow coefhcients will be zero here due to azimuthal symmetry. Clearly, this will no longer 
be true in the presence of azimuthal anisotropy in the detector acceptance. For example, let A{(j)) be the acceptance 
function characterizing the probability that a hadron is detected at azimuthal angle (0, with A{(j)) being normalized 
as [i,[ll 



A{(P)d(p = 2tt (16) 

In the lab fixed frame, the values of w„ in our method may be calculated from Eq.(3) (Sect. II). Note that for 
anisotropic detector acceptance there is an implicit factor of A{(j)) in the integrand in Eq.(3). One can then directly 
calculate the event average value v„ of these flow coefficients Vn, by averaging ?)„ over a large number of events (of 
similar type, e.g. centrality etc.). The actual value of Vn, for non-central collisions, is determined by the following 
expression 

C C A{4>)p^{c^)e'''"'^d4>)d^l; (17) 



27r . _^^2tt 



Here, ip denotes the orientation of the event plane in the lab fixed frame (see, Eq.(4) in Sect. II), which varies from 
event to event. We have made here explicit the acceptance function A{(j)) in the above expression. (For simplicity 
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we assume that A{(j)) does not depend on rapidity, pr etc.). For ideal detector with uniform acceptance A{(f)) = 1 
and by changing the order of d(j) and dip integrations one can easily see that u„ = as discussed in Sect. II (following 
Eq.(4)). However, for anisotropic acceptance, with A{(j)) ^ 1, Vn need not be zero. In fact, apart from a normalization 
(relating to angular average of p), w„ gives the Fourier coefficients for the series expansion of the function yl(0), thus 
allowing us to determine the acceptance function A{(^^ directly from experimental data [sl. ITsj. 

For the simple case when ^(0) is non-zero for every 0, one can directly use A(0), determined from Vn as explained 
above, to compensate for the effect of anisotropic detector acceptance by including additional factor of \j A{ip) in the 
integrand of Eq.(3) for calculation of v„. 



With this modification, Eq.(17) shows that the event average value S„ will again be zero, just as for the case of 
isotropic detector acceptance. The above Eq.(18) can then be used to determine the values of w™* as discussed in 
earlier sections. 

However, for more general case, e.g. for incomplete detector coverage, where ^(0) = for a range of the above 
simple method does not help. For such cases one can adopt the techniques discussed in [1, [l^ for the present case. 
Though, due to the important role played by the event plane determination, these techniques need to be suitably 
adopted for our case where a fixed lab frame is used for the analysis. In fact, it is more helpful to directly adopt 
the techniques used for CMBR case where one uses a fixed lab frame (say, galactic coordinates) for writing down 
temperature anisotropies in the sky in terms of spherical harmonics. For CMBR analysis also one needs to compensate 
for the effects of partial coverage of the sk y. T his happens either due to limited coverage of the sky by the detector, or 
due to galactic foreground as for WMAP [i6[. Adopting that approach for our case, we start with the experimentally 
determined flow coefficients v„ for each event (which, following [iq we call pseudo flow coefficients). 



1 

2^ 



■'^A(0)p^(0)d0 (19) 



where, now, A(<j)) may be zero for a range of ^ values. Using Fourier coefficients am for the acceptance function 
A{4)), we get. 



where vi is the true flow coefficient (defined with respect to the lab fixed frame). The above equation shows the 
mode- mode coupling of flow coefficients resulting from incomplete detector coverage [1, [l^ • For an approximate 
determination of vi we can use appropriate truncation of the above matrix equation to a finite set of linear equations 
and solve for w;. This analysis can be done for each event, and with the values of vi thus determined one can directly 
get the estimate of the true variance -D™'* for a set of events. 



VI. CONCLUSIONS 



In summary we emphasize the important lessons from CMBR analysis techniques. We have argued that important 
information about initial anisotropies of the system and their evolution in relativistic heavy-ion collisions can be 
obtained by plotting the root mean square values of the Fourier coefficients "D™* of the anisotropies in the fluctuations 
5p/p of the particle momenta, calculated in a fixed laboratory frame, starting from n = 1 upto large values of n ~ 30. 
(which corresponds to A ^ 1 fm). Note that n = 30 almost corresponds to wavelength of fluctuation A at the surface 
of the region, at r/,., being of order 1 fm. Fluctuations with wavelengths smaller than 1 fm presumably cannot 
be treated within hydrodynamical framework, so we restrict attention within this range of n. One will expect that 
beyond a critical value of n the nature of the curve should change in some qualitative manner indicating breakdown 
of underlying hydrodynamical description for smaller modes. The wavelength corresponding to that critical value 
of n will determine the smallest scale below which hydrodynamical description is not valid. The plot of -D™^ for n 
larger than this critical value will probe fluctuations in parton density at even smaller length scales and may provide 
a bridge with the perturbative regime. 

For non-central collisions our technique provides a direct method to probe elliptic flow by determining {i™* in lab 
fixed frame. If random fluctuations do not dominate then our results show that t)™" can directly probe |f„| for n ^ 2 



10 



also. This may be important for collisions of deformed nuclei. It is important to appreciate that a plot of v™^ tells us 
about the statistical nature of fluctuations and anisotropics of the initial plasma region (just as in the CMBR case), 
hence it will always have valuable information irrespective of its shape (e.g. a flat curve). Hydrodynamical simulations, 
with proper incorporation of such fluctuations can be used to predict directly the values of ■D™'* whose comparison 
with data will constrain/determine various physical inputs such as equation of state, viscosity etc. Further, in analogy 
with CMBR analysis, one can determine higher moments of momentum anisotropics (again in lab fixed frame, as we 
have discussed above) probing detailed nature of initial fluctuations, e.g. non-Gaussianity etc. We hope to present a 
more detailed analysis of such issues in a future work. We also intend to explore various aspects of horizon entering 
of density fluctuations which can be probed under laboratory conditions using RHICE. This may help in bringing 
at least some aspects of inflationary physics (e.g. causal aspects, excluding those features which relate directly to 
gravity) of the universe under some experimental control. 
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